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Abstract 

We use realistic pseudopotentials and a plane-wave basis to study 

the electronic structure of non-periodic, three-dimensional, 2000-atom 

(AlAs)„/(GaAs) m (001) superlattices, where the individual layer thicknesses 

n, m 6 {1, 2, 3} are randomly selected. We find that while the band gap of the 

equivalent (n = m = 2) ordered superlattice is indirect, random fluctuations 

in layer thicknesses lead to a direct gap in the planar Brillouin zone, strong 

wavefunction localization along the growth direction, short radiative lifetimes, 

and a significant band-gap reduction, in agreement with experiments on such 

intentionally grown disordered superlattices. 
PACS numbers: 73.20.Dx,78.66.-w,71.50.+t 
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Ordered semiconductor superlattices — produced routinely by epitaxial crystal growth 
techniques — are widely recognized for their unique electronic and optical properties To 
tailor the electronic properties (e.g., through band folding) one aims at growing ordered 
superlattices (o-SL) with definite values of the thicknesses n and m in (A) n / (G) m . One could, 
however, deliberately grow a disordered superlattice (d-SL) [0,0, where the individual layer 
thicknesses n,m,n' ,m' ,n" ,m" , . . . are selected at random according to a given probability 
distribution p a {n) (a = A, G). While the electronic structure of an o-SL is characterized by 
extended states and the formation of mini-bands, a truly one-dimensional disordered system 
(described, e.g., by the Anderson model) shows carrier localization and absence of dispersion 
gyj|. Sasaki et al. || grew ~1000 monolayer (ML) thick AlAs/GaAs d-SL's (i.e., A = AlAs, 
G = GaAs) with n and m chosen from a set of small integers {1, 2, 3}, viz., Pa = Pg = 
with p{l) = p{2) = p(3) = |. Since the average layer thickness is 2 monolayers, one can 
think of this d-SL as evolving from an (^4)2/ (GO2 o-SL by random substitution of A2 and G2 
layers by A\, A 3 , G±, and G 3 layers ( "5-doping" ) . Such disordered superlattices have shown 
surprising and unique optical properties relative to their parent o-SL [[$]: (a) strong and 
initially fast decaying (lifetime r = 0.25 ns at T = 77 K) photoluminescence (PL) intensity 
even though the equivalent o-SL has an indirect band gap and thus emits both weakly and 
slowly, (b) a large red shift (~60 meV) of the PL peak with respect to the equivalent o-SL, 
and (c) an order of magnitude slower rate of reduction of the PL intensity with temperature. 
These unusual properties of d-SL's appear very attractive for optoelectronic applications ||. 

In modelling the electronic structure of a d-SL []2]||[7|], one faces difficulties arising from 
the existence of two entirely different length scales: (i) The lack of translational symmetry 
requires the use of unit cells with a macroscopic length N m 1000 ML, equal to the total 
length of the d-SL (Nd ~ 300 nm, where d is the monolayer thickness), (ii) The spatial 
variations of the electron potential, however, occur on a microscopic length scale of about 
0.1 nm. While it is possible to rescale the microscopic length scale by replacing the periodic 
atomic potential by an external, rectangular potential this approach fails to describe 
the band structure (e.g., the indirect gap of the (AlAs)2/(GaAs)2 o-SL) in the present 
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regime of rapid layer fluctuations. To overcome the problems arising from the existence of 
two disparate length scales, we extended a microscopic pseudopotential description of the 
electron structure to a macroscopic length scale necessitated by the absence of translational 
symmetry. We use fixed (screened) atomic pseudopotentials that were carefully fitted to 
bulk band structures, effective masses, deformation potentials, band offsets, and energy 
levels in superlattices |§. The wavefunctions are expanded in about 30 plane waves per 
atom. For an N- monolayer superlattice along (001) with two atoms per monolayer the 
corresponding matrix dimension is therefore about 60iV x 60iV. Standard techniques to 
solve the Schrodinger equation require orthogonalization of the states of interest (i.e., near 
the band gap) to all lower-lying states. This leads to an N 3 scaling of the effort and becomes 
impractical for structures with iV rs 100 ML. We use instead the "folded-spectrum" method 
||, where eigenstates are obtained directly in an energy window of interest (e.g., near the 
band gap), without having to solve for any of the ~ 8N lower-lying eigenstates first, thus 
circumventing the need for orthogonalization. The effort scales linearly with N, allowing 
us to use the realistic, three-dimensional pseudopotentials, and to solve the Schrodinger 
equation in a highly flexible plane-wave basis even for N = 1000 ML. 

Application to d-SL's of AlAs/GaAs leads to an explanation of the large red shift, the 
enhanced oscillator strength, and the weak temperature dependence in terms of band-edge 
wavefunction localization. The source of localization is interesting: In truly one-dimensional 
disordered chains all states are in general localized |4],U . However, the laboratory-grown |3j 
d-SL's have extended layers in the (x, y) plane (perpendicular to the disorder direction z = 
[001]), so these quasi one-dimensional systems retain in fact a three-dimensional character, 
and the states need not be localized by disorder. We will show that the localization in the 
d-SL originates mostly from the formation of impurity-like, localized bound states due to 
insertion of 5-layers into the o-SL. 

In Fig. [L](a) we show the planar wavefunction average \i/je(z)\ 2 = J d 2 r±\ipE{f)\ 2 of a few 
occupied and unoccupied band-edge states of a d-SL with a 1000-monolayer unit cell. We 
see that these band-edge wavefunctions are sharply localized. We can quantify the effective 
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localization length (in monolayer units) for wavefunction ^ at energy E as || 

L e ff (^) = ^(El^l^)| 4 ) > (!) 

where the sum extends over the grid points i along z, and h is the number of grid points per 
monolayer. For a truly extended state, L e s is of the order of N, while for a state localized 
on one site, L eS = 1. We find that L eff < 15 ML for the band-edge states. The asymptotic 
decay length 7 _1 can be quantified by (\iJje(z)\) oc e~ 7Z , where the angular brackets denote 
averaging over the fast oscillations of ^^(z)! along z. Near the band edges, the calculated 
values are 7 « 0.2 ML -1 . 

We now investigate two possible mechanisms for the wavefunction localization appar- 
ent in Fig. |l](a): (i) chemical binding to 5-like "impurity layers" in an otherwise ordered 
(A) 2 /(G) 2 host, and (ii) a continuous increase of localization (measured by 7 or by L~^) 
with increasing degree of disorder. Scenario (i) is motivated by the fact that in one di- 
mension an attractive 5-potential always has one bound state, whereas scenario (ii) is valid 
for the Anderson model with on-site disorder obeying a continuous probability distribution 
0. In our case of a discrete probability distribution p(n), we connect the two pictures by 
starting with the reference o-SL (A) 2 /(G)2, and gradually substituting A 2 and G 2 layers by 
"wrong" layers with n = 1 or n = 3. We measure the degree of disorder by counting the 
relative frequency R of these layers, i.e., 

R = p(l)/p(2)=p(3)/p(2). (2) 

The fully disordered SL has R = 1, the single 5-layer in an ordered (A) 2 /(G) 2 superlattice 
corresponds to R ~ N~ x , while the perfect o-SL has R = 0. 

To understand the possibility of impurity-like localization, consider for example a G3 5- 
layer embedded in the otherwise perfect o-SL ■ • ■ A2G2A2G2A2G2 ■ ■ ■, thus converting it into 
• ■ • A2G2A2G3A2G2 ■ ■ ■■ If the G3 5-layer is attractive to electrons (holes) it will bind a state 
below the conduction-band minimum (CBM) [above the valence-band maximum (VBM)] of 
the o-SL We find that a (GaAs)3 5-layer in the (AlAs)2/(GaAs)2 o-SL indeed binds 



an electron and a hole [Fig. |](b)], while an (AlAs)3 layer binds an electron but does not 
bind a hole [Fig. [j](c)]. The bound-state localization lengths 7 _1 and L c r obtained with a 
single 5-layer are similar to those obtained in the fully disordered SL (Fig. [I]), suggesting 
that the same mechanism of localization could be at work in both cases. As one increases 
the concentration R of randomly positioned 5-layers, one finds more bound states which are 



eventually forming a quasi continuum inside the band gap jy] . This is illustrated in Fig. y, 



where the band-edge energies of d-SL's with N = 128 ML [12[] are plotted as a function of the 
degree of disorder R. As R approaches zero, the band edges merge with the finite binding 
energies of an isolated 5-layer [scenario (i)], and not with the unperturbed band edges of the 
o-SL [scenario (ii)]. This significant observation suggests that the localization energy in the 
d-SL comes mostly from impurity binding. Figure |2] shows that for large degrees of disorder 
R, the band edges are pushed further into the gap. To isolate the effect of pure disorder 
from the effect of impurity-like localization, we also show in Fig. |2| the band edges of an 
o-SL containing a periodic array of 5-layers of concentration R, i.e., separated by a distance 
~ Br 1 . We see that even for an array of closely spaced 5-layers (R — > 1), the binding energy 
does not increase, indicating a negligible interaction between the neighboring, coherently 
arranged bound states. In contrast, in a d-SL the 5-layer-like bound states are arranged 
incoherently , leading to a band tail. These studies show that the localization length of the 
band-edge states in a d-SL is decided already by the chemical, impurity-like binding of a 
single 5-layer [Figs. |](b) and [l](c)] . Furthermore, the energy position of the gap levels at 
small to intermediate degrees of disorder R is also determined by the properties of non- 
interacting, periodically arranged 5-layers (Fig. |2|). For higher values of R, disorder shifts 
and broadens the gap levels further into the gap, without modifying their localization length. 

To illustrate the dependence of the gap-level shifts on the particular form of disorder, we 
compare the results obtained for the d-SL with those found in a partially-ordered superlattice 
(po-SL): Arent et al. [|13[ grew po-SL's with the same distribution function p(n) for the A n 
layers as in a d-SL, but preserved long-range order by requiring that each A n layer be 
followed by a G^-n layer. Therefore, at positions 1,5,9,. . .there is always an A layer, and at 



positions 4,8,12,. . .there is always a G layer, while at the "sandwiched" positions, A and G 
layers are distributed randomly. We see in Fig. |2] that the presence of long-range order leads 
to an even larger shift of gap levels than in the equivalent d-SL at the same R value. The 
band-edge wavefunctions have the same characteristic L e ff and 7 as in a d-SL. Thus, absence 
of long-range order (in the d-SL) is not essential for obtaining large band-edge shifts. 

The calculated band-gap reduction of the d-SL and po-SL with respect to the o-SL 
(Fig. |]) is consistent with experiment p|Jl3||: We find band gaps of 2.09, 1.94, and 1.87 eV 
for the o-SL, d-SL, and po-SL, respectively, compared with the experimental PL emission 
lines at 2.02, 1.96 0, and 1.87 eV @, respectively. 



To determine if the d-SL has a direct or indirect gap in the two-dimensional Brillouin 
zone, we show in Figure |3| the dispersion of the band-edge states of the d-SL (solid lines), 
o-SL (thin lines) and single 5-layer (dotted lines) along the symmetry lines £ and A, i.e., 
from f to M = ^(1, 1) and from F to X = ^(1, 0), respectively. The thin horizontal lines 
denote the band edges of the underlying o-SL. We find that the conduction bands of the d-SL 
dip below these lines. The difference ("binding energy") increases in the order M — > F — > X. 
|T4| Thus, the large binding energy at T pulls the conduction-band edge below the one at 



M by 60 meV, making the d-SL a direct-gap material ||15||, even though the o-SL is indirect 



(with CBM at M). This suggests strongly that the observed || strong PL intensity is due 
to the occurrence of a direct band gap, leading to efficient recombination of electrons and 
holes localized in the same spatial region along z [see, e.g., the same positions along the 
chain of states v2 and c3 in Fig. 0(a)]. The localization along z relaxes the fey-selection rule, 
thus further enhancing the oscillator strength. The enhanced oscillator strength is reflected 
by short radiative lifetimes r: we calculate r = 1 ns for the v2— >c3 transition at energy 1.96 
eV [Fig. [l|(a)]. These radiative lifetimes are 1000 x faster compared to those measured in 
indirect-gap o-SL's (r « 5.5 /is at T = 2 K) JT§. 

To understand why the PL intensity in a d-SL has a weaker temperature dependence 
than in an o-SL, consider vertical transport to non-radiative recombination centers. This 
channel will be inhibited, unless there are strongly overlapping electron states within an 
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energy ~ kT of each other. Strong overlap occurs when ((E,T) = /cTp(.&r)L e ff(.&r) > 1, 
where p{E) is the one-dimensional density of states, and is a typical energy that is 
thermally occupied at temperature T. Figure || shows p(E), L e s(E), and ((E,T) obtained 
by averaging over 100 realizations of N = 2000 ML d-SL's, calculated within the Kronig- 
Penney approximation with band offset fitted to our pseudopotential calculations. We find 
that: (i) the DOS exhibits a ~200 meV wide band tail extending below the CBM of the 
o-SL (the zero of energy); (ii) L e g is almost constant (~20 ML) in the band tail; (iii) at the 
thermally populated levels [denoted by the arrows in Fig. f|(b)] C ^ 1> thus vertical hopping 
is suppressed. Consequently, the d-SL will have a weaker temperature dependence of the 
PL decay, because higher temperatures will be needed to dissociate the electron-hole pairs 
in the vertical dimension. 
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FIGURES 

FIG. 1. Planar average of the wavefunctions squared in (a) a 1000-monolayer d-SL, (b) an 
o-SL host containing a single (GaAs)3 5-impurity and (c) a single (AlAs)3 <5-impurity. Unoccupied 
states are labeled cl, c2,. . . ; occupied states are labeled vl, v2,. . . , and are plotted in the negative 
direction with a small offset for clarity. L e g and 7 are given in ML and ML -1 , respectively. 

FIG. 2. Band-edge energies of disordered (d-SL, diamonds), partially ordered superlattices 
(po-SL, pluses), and a periodic array of 5-layers embedded in an o-SL (thin horizontal lines) as 
a function of R [Eq. @]. Lines are guides to the eye. The thick horizontal lines denote the 
unperturbed band edges of the parent o-SL. The vertical bars on the d-SL data points denote the 
range of binding energies obtained from ~10 different realizations of a d-SL with length N = 128. 

FIG. 3. Dispersion of the band-edge state in a d-SL (solid lines), o-SL (thin lines), and a single 
5-layer in an o-SL host (dotted lines) along the A and S symmetry lines in the planar Brillouin 
zone. Also shown are the (OOl)-projected bands of the o-SL (shaded area). The supercell size was 
N = 128. 

FIG. 4. (a) Density of states p{E) (solid line), L g r{E) (dashed line) averaged over 100 realiza- 
tions of N = 2000 ML d-SL's, calculated with the Kronig-Penney model. Dots indicate L e f[(E) ob- 
tained for an iV = 128 ML pseudopotential calculation, (b) The product (,{E, T) = kT ' p{E)L c ^{E) . 
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Figure 4 (revised), Mader and Zunger 
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